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A very small dispersion in the speed of light may be observable in Fermi time- and energy-tagged 
data on variable sources, such as gamma-ray bursts (GRB) and active galactic nuclei (AGN). We 
describe a method to compute the size of this effect by applying the Feynman sum-over-histories 
formalism for relativistic quantum electrodynamics to a discrete model of space-time called causal 
set theory. 



I. INTRODUCTION 

The vacuum speed of light may be a function of 
photon energy. Even a very small such dependence 
might yield energy-dependent time delays which, over 
the long journeys from distant astronomical sources, 
may accumulate to a level measurable by the Fermi 
Gamma Ray Space Telescope [ ]. The expectation 
from general considerations [2] is that the scale of the 
effect is linear in the photon's energy: 



dt dE 

T Epi anc k 



where the Planck energy is 



E 



Pla 




1.22 x 10 19 GeV . 



(1) 



(2) 



This relation gives delays over astronomical distances 
of seconds (GeV energies) to hours (TeV energies). 

A motivation for this suggestion is the possibility 
that space-time may be lumpy, fuzzy, or even discrete, 
on small scales, and that this lumpiness may affect 
photon propagation in a dispersive way. We adopt 
the extreme viewpoint provided by causal set theory, 
an approach to quantum gravity pioneered by Rafael 
Sorkin [ ], with the founding principle that space- 
time is a partially ordered set of purely discrete points. 
The partial order nicely reflects the causal relation- 
ships between events fundamental to special relativ- 
ity. We also adopt the Feynman sum-over-histories 
quantum electrodynamics, which is greatly simplified 
by the postulated discreteness: path integrals, and 
the difficulties associated with divergences and defin- 
ing the measure on the space of all paths are replaced 
by discrete sums. 



II. CAUSAL SETS 

Causal set (CS) theory postulates that at very small 
scales - below those of ordinary experience, and even 



below the smallest scales accessible to particle physics 
experiments - space-time consists of a set of discrete 
points. The usual space-time continuum is thus re- 
garded as a macroscopic construct that does not exist. 
One only requires that the large scale limit of the CS 
recovers known, macroscopic, continuum physics. 

In mathematical terms, a causal set is a pair (C, -<), 
where C is a set and ^< is a binary relation on C satis- 
fying the following properties (x, y, z etc. are general 
points in C): 

1. Partial Order: Some, but not all, pairs x,y 
of points are ordered: x -< y\ 

2. Transitivity: If x -< y and y -< z, then x -< z; 

3. Anti-reflexivity: No point is related to itself: 

4. Local finiteness: For all x, y E C, the set 

[x, y] = {z E C: x -< z and z -< y} is finite. 

The first two items are consistent with, and are taken 
to express, the causal relations of special relativity: 
x < y (read "x precedes y") means that y is in the 
forward light-cone of x, and can affect x. If sepa- 
rated by a space-like interval x and y are not or- 
dered. Transitivity here simply expresses the relation 
between nested light cones. 

One consequence of the above ansatz is that all 
space-time information is contained in the causal con- 
nections among the points. For example, any topo- 
logical or metric notions must be extracted from 
the discrete (i.e., combinatorial) structure of (C, -<) 
[6, 7, 10, 11], without any reference to a continuum. 
This also means that the action (see below) for the 
step x — > y can depend on only the causal relation 
between the two points, plus the energy and polariza- 
tion of the photon. We make use of this simplifying 
fact below. 

Note that arbitrarily defined partial orderings 
among the points in general can yield causal sets that 
cannot be realized as a subset of any Lorentzian man- 
ifold. The sprinkling procedure [5] for generating CS 
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order relations was invented to avoid these micro- 
scopic inconsistencies, to avoid other problems that 
might foul up the macroscopic continuum limit [14] 
mentioned above, and to ensure a kind of statistical 
Lorentz invariance. The process begins by randomly 
sprinkling points in a continuous Lorentzian space- 
time manifold M, of the appropriate dimension. For 
example, to model flat space, the points would be uni- 
formly and independently distributed[15] in M. The 
fundamental density g of this distribution is frequently 
taken to be such that the mean (continuum) separa- 
tion of points is the Planck scale. Then the partial 
order relations among the sprinkled points are derived 
from the causal (light-cone) relations assessed in the 
continuum M. Here we adopt this sprinkling concept 
as defining causal sets, which in the literature are said 
to be embedded into the continuum M. 



III. PHOTON PROPAGATION: FEYNMAN 
SUMS OVER SPACE-TIME PATHS 

Given a causal set (C, -<) embedded in some M — 
R d+1 , where d is the spatial dimension and there is one 
time dimension, we consider photon propagation in C. 
Given spacetime points x,y, one asks for the value of 
the propagator K(x,y), i.e. the quantum- mechanical 
amplitude for a photon emitted at x to be detected 
at y\ the probability of this transtion is \K(x,y)\ 2 . 
The standard sum over paths formalism [3, 4] states 
that the amplitude for each path from x to y is the 
product of the amplitudes of each step comprising the 
path, and the amplitude for the process is obtained 
by summing these amplitudes over all paths from x 
to y. 

It is more realistic to take <S, a subset of C, to repre- 
sent an unresolved astronomical source, and another 
subset V to represent the detector. Then the usual 
prescription is that the amplitude for the process of 
interest - photons emitted from the source arriving at 
the detector, symbolically V - is the sum of the 
amplitudes over all paths starting beginning at some 
point in S and ending at some point in V. The proba- 
bility of the process V is the square of the complex 
amplitude of the process. This prescription is valid if 
and only if one is not able to determine which point in 
V the photon arrives at, nor from which point in S the 
photon originates. If the detector were to be divided 
into pixels, such that the identity of the pixel at which 
the photon arrives is knowable, then as usual ampli- 
tudes are additive within pixels, but probabilities are 
additive across pixels. 

To compute the propagator K(x, y) we use the Feyn- 
man path integral approach which asserts that the 
propagator is the "integral" of expos' (7) /ft) over the 



space of all continuous paths 7 from xtoy, where 5(7) 
is the action associated with 7. Whereas in the contin- 
uous case this integral is not well-defined in a rigorous 
sense and its meaning, definition and computation are 
still an area of investigation, in the causal set scenario 
the integral reduces to a sum, to which it is easier to 
give a precise meaning. Related work can be found in 
[8, 9], but unlike these references our approach takes 
the all paths spirit of [4] quite literally, allowing causal, 
non- causal, and superluminal paths and relying on the 
phase averages for large action paths as appropriate 
weights. 

All that remains is specification of the action S(x, y) 
for an elementary step x — > y, for arbitrary points x 
and y in C, for the corresponding amplitude is [ ] 

S(x,y) = S e iii ^ 1 (3) 

where h is Planck's constant, and So is a normalization 
factor which can be determined by considering the 
operator for an infinitesimal advance in time ([3], §6). 

IV. MATHEMATICAL RESULTS 

To effect a suitable choice of the action just de- 
scribed, take the amplitude of each causal or non- 
causal jump to be a or 6, respectively, where a and 
b are suitably chosen complex numbers. Given a path 
7 = X0X1 - — x n consisting of n steps, set £(7) = 
a n ~ k b k , where k is the number of non-causal steps 
in 7, and n — k the number of causal ones. The prop- 
agator from x = xo to y = x n is then 

00 n 

K (*> f) = EE a n - k b k Pn , k (x, y), (4) 

n=l A;=0 

where p n ^(x : y) denotes the number of paths of length 
n from x to y with exactly k non-causal steps. We 
now adopt a statistical view and, in a slight abuse 
of notation, treat K and p n ,k( x ,y) as averages over 
random sprinklings. 

One obvious problem with this definition is that for 
n > 2 and k > 1, Pn,k(%,y) is infinite. To rectify this 
problem, we pick a sufficiently large number r > and 
allow jumps (both casual and non-causal) only up to 
distance r. That is, 7 = xoX\ • • • x n is defined to be an 
admissible path only if each step satisfies d(xj, Xj+i) < 
r, where d is the distance metric. This distance should 
be defined in the CS sense (e.g. as in [ ]), but in this 
preliminary study we have ignored this nuance. If we 
denote the corresponding propagator by K r , then 

00 n 

K r (x, y) = J2H a n ~ k b k p r n, k ( x > V)> ( 5 ) 

n=l k=0 
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where p r n k (x, y) is the expected number of admissible 
paths of length n from xtoy with exactly k non-causal 
steps [we will call such paths (n, fc)-paths]. One can 
think of the sum in equation (5) either as an ordinary 
series (whose convergence will be assured if \a\ and 
\b\ are sufficiently small) or as an oscillatory series, 
whose convergence is defined differently. We choose 
the former interpretation of equation (5). 

To compute K r (x,y) for any x,y G R d+1 , we need 
a suitable expression for p r nk (x,y). The amplitudes 
of causal and non-causal steps will be represented by 
defining two functions v r ,jJL r : R d+1 — > R as follows: 



v r {z) 



1 if z < and d(z, 0) < r 
otherwise, 



(6) 



1 l ± are the characteristic functions of the future and 
past light cones of the origin. See Figure 2. 




and 



fJL r (z) 



1 if z ^ and d(z, 0) < r 
otherwise. 



(7) 



See Figure 1. Observe that v r {x — y) = 1 if and only 
if there is an admissible non-causal jump (i.e., path 
of length one) from x to y, and \i r (x — y) = 1 if and 
only if there is an admissible non-causal jump from 
x to y. It is not hard to see that v r = Ib^L- and 
H r = 1 B r 1l + where 1 B r is the characteristic function 
of the ball B r of radius r centered at the origin, and 



FIG. 1: The region marked v r is the backward light 
cone of 0, truncated at distance r. The remaining region 
marked fi r includes the forward light cone and the remain- 
ing region with space-like separations from, and therefore 
causally disconnected from, 0. 

It is possible to derive explicit formulas for 
Pn k( x i y) usm g convolution-type integrals of functions 
v r ,ii r . For lack of space, we do not present these for- 
mulas. It can then be shown that the functions p r n k 
satisfy the following family of integral equations: 



Pn+l,k( X > V) = Q I Pn,k( X ' Z )"r(z ~ V) dz + Q / P^k-lfa z)fl r (z ~ y) dz, (8) 

JRd+l jRd + 1 

for all n > 1 and 1 < k < n. A lengthy calculation yields the following integral equation for the propagator: 

ag / K r (x, z)u r (z — y) dz + bg / K r (x, z)fi r (z — y) dz — K r (x, y) — av r (x — y) — b/j Jr (x — y) (9) 

CO CO » 

- Y) bn Pn,n( X i y)+ b QyZ bn / Pn,n( X > Z )^r(z - y) dz. 



n=2 



n=l 



To solve equation (9), we observe that both K r (x, y) 
and p^ n (x,y) are translation invariant. This means 
that Kr(x,y) = ^ r (x - y) and p r nn (x,y) = P r n {x - 
y), for some functions ^ r ,^n ^ ^ 1 (R d+1 ). Note also 
that the integrals in equation (9) are convolutions. 
Applying the Fourier transform^, solving for \I/ r , and 
applying the inverse Fourier transform yields: 



bjJL r ) 



*(n+l) 



(10) 



n=0 



lution power of av r + b\i r . To compute the geometric 
convolution sum in equation (10), we need to recall 
that the space L 1 (R d+1 ) of Lebesgue integrable func- 
tions is a Banach algebra [ ] with respect to multi- 
plication given by convolution *, but that it does not 
possess a unit element. However, a unity - usually 
denoted by 5 (one can think of it as the Dirac delta 
"function") - can be adjoined to L x (R d+1 ) to obtain 
a new Banach algebra L x (R d+1 ) with unity (cf., [12], 
§10.1). Then the sum in equation (10) becomes 



where (av r + 6/i r )*( n+1 ) denotes the (n + l)-st convo- 



\I/ r = (av r + bfi r ) * {5 — g(av r + bfi r )} x , (11) 
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FIG. 2: Jumps from x to y and x to z are admissible, but 
not from x to w. 



where the inverse is taken in the extended algebra 

Now let S and V be subsets of R d+1 representing 
the source and the detector (see Section III) . Then the 
probability a photon traveling from S to V is 



Prob(<S V) 



l r (x, y) dxdy 



(12) 



We leave the discussion of this integral for various S 
and V to our forthcoming longer paper. 



V. NUMERICAL RESULTS 

A matrix formalism simplifies numerical evaluation 
of the propagators for small CSs. For a causal set 
consisting of N points x n , n = 1, 2, ... N, define the 
amplitude matrix 



(13) 



generalizing the connection matrix of graph theory by 
replacing the 1 — indicator of connectivity with the 
full amplitude from eq. (3) for the step between the 
points. The product of A with itself, 



gives the amplitude for two-step transitions n — > m 
with a single intermediate stop at j. Similarly [A fc ] n?m 
gives the amplitude of the same transition with k — 1 
intermediate stops, and the series 



(15) 



gives the amplitude summed over all paths. As long 
as || A || < 1 for some norm, this series converges to the 
computationally useful form 



Q = (I-A)" 



I 



(16) 



2 1 - V 4 A- 



where I is the indentity matrix. 

So the amplitude of the transition S — > V is the 
sum of the elements of Q in the sub-block with column 
indices labeling points in S and row indices labeling 
points in D, and the final probability for S — > V is the 
absolute square of this sum. 

Considering this result and the comments about 
pixels which identify detected photons made above, we 
see how precisely and simply these matrix operations 
implement the fundamental sum-over-histories quan- 
tum rules: multiply amplitudes of steps on a path; 
add amplitudes of all paths between indistinguishable 
starting and final states; add probabilities correspond- 
ing to distinguishable final states. 

After this construct was developed during a CAM- 
COS project at San Jose State University we discov- 
ered that Steven Johnston has a similar approach [8]. 

The formalism developed here seems to be a promis- 
ing approach to computing photon dispersion in the 
context of causal set theory. A future publication will 
elaborate the analytical approach in Section IV and 
the numerical approach in Section V, including results 
on the speed of light as a function of energy. 
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